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Abstract. Taking the total angular momentum conser- 
vation in place of the tangential momentum conservation, 
and considering the square and higher power terms of or- 
bital eccentricity e, the changes of orbital elements of bi- 
naries are recalculated for wind accretion scenario. These 
new equations of orbital elements are used to calculate the 
properties of barium stars. Our results show that, during 
the evolution of a binary system, the system widens as it 
loses mass, and the orbital period increases, while orbital 
eccentricity remains nearly constant, which can quantita- 
tively explain the observed (e,logP) properties of normal 
G, K giants and those of barium stars. The results reflect 
the evolution from G, K giant binaries to barium bina- 
ries, namely, the orbits of barium stars have been modified 
by the mass-transfer process responsible for their chemi- 
cal peculiarities, whereas most of the G, K giant binaries 
are probably pre-mass transfer binaries. Moreover, the re- 
sults showed that the barium stars with longer orbital 
period P>1600 days may be formed by accreting part of 
the ejecta from the intrinsic AGB stars through wind ac- 
cretion, while those with shorter orbital period may be 
formed through dynamically stable late case C mass trans- 
fer or common envelope ejection. 

Key words: Stars: late-type- Stars: chemically peculiar- 
Stars: AGB- Stars: Mass loss- Binaries: spectroscopic 



1. Introduction 

Extrinsic AGB stars include various classes of G- and 
K-type barium stars and those cooler S and C stars in 
which ^^Tc (ti=2x10^ yr) is not observed. It is gen- 
erally believed that they belong to binary systems and 
their heavy-element overabundance come from accretion 
of the mass ejected by their companions (the former AGB 
stars, now evolved into white dwarfs) (McClure et al. 1980; 
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Boffin & Jorissen 1988; Jorissen & Mayor 1992; Jorissen 
& Van Eck 2000; Jorissen et al. 1998; Bohm-Vitense et 
al. 2000; Karakas et al. 2000). The mass exchange took 
place about 1 x 10^ years ago, so the ^^Tc produced in 
the original thermal pulse AGB (hereafter TP- AGB) stars 
have practically all decayed. The accretion may either be 
disk accretion (Iben & Tutukov 1985) or common envelope 
ejection (Paczynski 1976). Han et al. (1995) detailedly in- 
vestigated the evolutionary channels for the formation of 
barium stars. 

Boffin & Jorissen (1988) calculated quantitatively the 
variations of orbital elements caused by wind accretion in 
binary systems. They also estimated the heavy-element 
overabundance of barium stars. In a subsequent paper 
(Boffin & Zacs 1994), the heavy-element overabundance 
of barium stars was calculated using the same technique, 
and the correlation with orbital period was interpreted. 
The calculated results were compared with the observa- 
tions (Zacs 1994). In our calculation of the overabundance 
(Chang et al. 1997), we considered the variation in the 
binary separation, dr ^ 0, and combined the nucleosyn- 
thesis scenario of intrinsic TP-AGB stars (Zhang et al. 
1998) with the scenario of successive pulsed accretion and 
mixing of wind mass. 

Some important conclusions have been obtained in the 
theory of wind accretion, but previous calculations on or- 
bital elements were not very reliable because of the ne- 
glect of the 6r/r term and the adoption of tangential mo- 
mentum conservation (Chang et al. 1997 and references 
therein). For the rotating binary system with wind mass 
loss, total angular momentum conservation is a more real- 
istic assumption than tangential momentum conservation. 
Also, for the sake of simplicity, previous calculations re- 
garded orbital eccentricity as e ^ 1 (Huang 1956; Boffin 
& Jorissen 1988; Theuns et al. 1996). Observations show 
that the eccentricities of almost half of all barium stars 
are too large to be neglected, with the largest observed 
value = 0.97 (Jorissen et al. 1998). Thus the assumption 
of e <C 1 for the barium star systems is not valid. 
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The (e,logP) diagram is a very useful tool to study bi- 
nary evolution given the abundant information contained 
in it. For example, the distributions of orbital period P 
with the eccentricity of progenitor systems correlate with 
the final characteristics of binary systems through mass 
accretion. Generally, barium stars evolved from normal 
red giant binary systems. Observations show that the or- 
bital eccentricities of Ba stars are smaller than those of 
the G, K giants from open clusters, and the (e,logP) dia- 
grams of the two systems show strong correlation (Joris- 
sen & Boffin 1992; Boffin et al. 1993; Jorissen et al. 1998; 
Jorissen 1999). Few quantitative calculations have been 
put into this correlation though they are very important 
to understanding the binary evolution. 

In this paper, we adopt the hypothesis that the total 
angular momentum is conserved, and we do not neglect 
the square and the higher power terms of eccentricity, in 
order to recalculate the change equations of orbital el- 
ements of binary systems caused by stellar wind accre- 
tion. Then these new equations are used to calculate the 
properties of barium stars produced by this mechanism. 
Hence we quantitatively explain the observed relations be- 
tween the (e,logF) diagram of normal G, K giants and 
that of barium stars. In Sect. 2, we analyze the proper- 
ties of (e,logP) diagrams of barium stars and normal G, 
K giants. In Sect. 3, we present the angular momentum 
conservation model of the wind accretion. Our results and 
analysis are illustrated in Sect. 4. And in Sect. 5, we draw 
some conclusions and discuss the results. 

2. The (e,logP) diagram 

The observations of the G, K giants from open cluster and 
barium stars arc the same as in Jorissen ct al. (1998), and 
are given in our Fig. 1 and 2 respectively. Analyzing the 
observations, we can understand the properties of them as 
following. 

(1) At any given orbital period, the maximum eccen- 
tricity found among barium systems is much smaller than 
for cluster giants. And the average eccentricity of barium 
stars is generally smaller than that of cluster giants. The 
average orbital eccentricity of barium stars is 0.14, while 
the corresponding mean value of cluster giant is 0.23 (Bof- 
fin et al. 1993). 

(2) In the two (e,logP) diagrams, the maximum ec- 
centricity line of barium system is almost parallel to the 
corresponding maximum eccentricity line of cluster giants. 
It seems that the maximum eccentricity line of cluster gi- 
ants is moved in a parallel way with the increase of orbital 
period to reach the barium systems. Boffin et al. (1993) 
exhibited this character. 

(3) In the (e,logP) diagrams, the average eccentricity 
line of barium stars is almost parallel to the corresponding 
line of cluster giants. It seems that the latter is moved in 
a parallel way to the longer period region to reach the 
barium systems. 



(4) In the (e,logP) diagram of normal G, K giants, 
the maximum orbital period for a circular orbit is ~350 
days; while the maximum orbital period for a circular orbit 
is ^2000 days in the (e,logP) diagram of barium stars, 
which is longer than the corresponding threshold of cluster 
giants. The (e,logP) diagrams of S and CH stars are very 
similar to those of barium stars (Jorissen et al. 1998). 

(5) Although, at a given orbital period, the maximum 
eccentricity found among barium systems is much smaller 
than for cluster giants. Still, the (e,logP) diagram of bar- 
ium systems shows that quite large eccentricities are found 
among barium stars, and the highest value is up to 0.97 
(HD 123949), yet at large periods (P~9200 days for HD 
123949). 

The reasonable binary evolutionary model of orbital el- 
ements should quantitatively explain the above mentioned 
observational properties. 

3. The angular momentum conservation model of 
wind accretion scenario 

For the binary system, the two components (an intrinsic 
AGB star with mass Mi, the present white dwarf star, and 
a main sequence star with mass M2, the present barium 
star) rotating around the mass core C, so the total an- 
gular momentum is conserved in the mass core reference 
frame. If the two components exchange material through 
wind accretion, the angular momentum conservation of 
the total system is showed by: 

Aii^r'^e) = ujrfiAMi + AM2) + r2v{AMi + AM2), (1) 

where fj, is reduced mass, and r is the distance from M2 to 
Ml . ri , r2 are the distances from Mi , M2 to the mass core 
C respectively. a;(=27r/P, P=2ttA'^{1 - e'^)i /h, is orbital 
period) is angular vclocdty. v is an additional effective ve- 
locity defined through the angular momentum variation 
in the direction of orbital motion of component 2. The 
first term on the right side of the equal-sign is the an- 
gular momentum lost by the escaping material, and the 
second term is the additional angular momentum lost by 
the escaping material. 

Using the methods similar to those adopted by Huang 
(1956) and Boffin & Jorissen (1988), but considering the 
angular momentum conservation of the total system and 
not neglecting the square and higher power terms of ec- 
centricity, we can obtain the formulas for the changes in 
the orbital elements: 
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where ^ is the semi-major axis of the relative orbit of 

component 2 around 1, and e is the orbital eccentricity 
(more details can be found in Appendix). Here, we take 
t; = (Boffin & Jorissen 1988). 

For the mass being accreted by the barium star, we 
use the Bondi-Hoyle accretion rate (Theuns et al. 1996; 
Jorissen et al. 1998): 



a 



GM2 



1 



1 + {Vovh/Vcjy 



AMi, 



(4) 



where a is a constant expressing the accretion efficiency, 
and a~0.05 in the situation of interest here according 
to the detailed hydrodynamic simulations (Theuns et al. 
1996; Jorissen et al. 1998). e is the orbital eccentricity. Vej 
is the wind velocity and Vorh is the orbital mean veloc- 
ity. After fixing the initial conditions, for the mass AMi, 
ejected at each pulse by the primary star, we can solve 
the Eqs. (2) — (4) for AM2, the mass accreted by the sec- 
ondary star. The final mass and critical Roche lobe radii 
of the intrinsic AGB stars, and the initial envelope mass of 
the extrinsic AGB stars are all calculated using the same 
formulas as in Boffin & Jorissen (1988) in which the mass 
ratio is q = M1/M2. 

4. Results cind cinalysis 

Adopting the wind accretion scenario above mentioned, 
we calculate the orbital elements of barium systems. 
We take as standard case: Mi^o=3.0 Mq, M2fl=l.2>MQ, 
i)ej=15 km s~-^ (Boffin & Zacs 1994). The results are 
showed in Figs. 1—4 and Table 1. 

According to the suggestions of Jorissen et al. (1998) 
and Zhang et al. (1999), the barium stars with longer 
orbital periods (P>1500 d and P>1600 d respectively) 
are formed through stellar wind accretion scenario. So we 
adopt 1600 days to be the lower bound of wind accretion 
for the formation of barium stars in our calculation. Fig. 1 
show the maximum eccentricities line of the normal giants 
(solid line), which is regarded as the initial orbit of the bi- 
nary system. The observations are for the G, K giants in 
open cluster (Mermilliod 1996; Jorissen et al. 1998). While 
the solid line in Fig. 2 represents the calculated eccen- 
tricity of barium stars through our wind accretion model, 
which is near to the observed maximum eccentricities line 
of barium stars. 

Fig. 3 show the (e,logP) results about mean orbital ec- 
centricities. The dashed line (1) represents the fit of mean 
eccentricity of observed cluster giants, which is taken from 
Boffin et al. (1993); the dashed line (2) represents the fit 
of mean eccentricity of observed barium stars taken from 



(3) 



1 






0.9 


G-KIII 




0.8 










/ • • 






/• . • 


0.7 










• / 


0.6 


• 




05 0.5 


• 


/ ••. • • 


0.4 


• 








0.3 






0.2 










0.1 








• ••• , 





r ' •mwmk— * « 





100 



1000 

p 



10000 



Fig. 1. The (e,logP) diagram for the G and K giants. The 
solid represents the maximum orbital eccentricity line of 
red giant binaries. 
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Fig. 2. The (e,logP) diagram for barium stars. The solid 
represents the calculated maximum orbital eccentricity of 
barium star systems created by the wind accretion sce- 
nario. 



Boffin et al. (1993). The thin solid line (4) represents the 
orbital eccentricity of barium stars with 77 = 0.685 in Eq. 
(11) of Bofiin et al. (1993) (Bofiin & Jorissen (1988) pre- 
dicted values of -q in the range 0.3 to about 1). In calcula- 
tion, we firstly fix the eccentricity of barium stars at given 
periods (thin solid line (4)), then use our wind accretion 
model to determine the corresponding initial eccentricity 
of binary systems (thick solid line (3)), which nearly over- 
lap the observed mean values of cluster giants (dashed line 
(1))- 
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Fig. 3. Mean eccentricity— logP diagram for binary sys- 
tems. The dashed line (1) represents the fit of mean eccen- 
tricity of observed cluster giants taken from Boffin et al. 
(1993), and the dashed hne (2) represents the correspond- 
ing fit of mean eccentricity of observed Ba stars. The thin 
solid line (4) represents the orbital eccentricity of barium 
stars with rj = 0.685 in Eq. (11) of Boffin et al. (1993), 
and the thick solid line (3) represents the corresponding 
initial eccentricity of binary systems. 

Table 1. Variations of the eccentricity e during the process 

of mass loss by wind as a function of the mass ratio q for the 
different initial eccentricities. The initial eccentricities are 0.1, 
0.3, 0.5 and 0.75 respectively. 
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Fig. 4. Evolution of the Roche lobe radius Rr in AU 
(solid lines) for the different final orbital periods during 
the process of mass loss by wind as a function of the mass 
ratio q. The final values of orbital periods (in days) are 
indicated. 

Those results for the (e,logP) diagrams show that the 
barium stars indeed evolve from the G, K giants, and 
the properties of barium stars with longer orbital peri- 
ods (P>1600 d) can be explained by the wind accretion 
scenario. 

Table 1 reports the variations of the eccentricities e 
during the process of mass loss by wind as a function of 
mass ratio q {q=Mi/M2) for different initial eccentricities 
of binary systems. The initial eccentricities are 0.1, 0.3, 
0.5 and 0.75 respectively. The results show that e is very 
nearly constant in the course of the wind accretion, which 
is different from the result of Theuns et al. (1996). 

Fig. 4 shows the evolution of the Roche radius due 
to mass loss via the wind as a function of the mass ra- 
tio q=Mi/M2 for different final orbital periods (in unit 
of days) of barium systems. The results show that, for a 
given final period of a barium star, the mass loss through 
stellar wind leads to the increase of the Roche radius, thus 
making Roche lobe overflow (hereafter RLOF) more diffi- 
cult to achieve. Moreover, the longer final orbital period, 
the larger Roche radius will be obtained, hence it is more 
difficult to occur the Roche lobe overflow scenario, so that 
the circular orbit forms with more difficulties. 

5. Conclusions cind discussion 

Taking the conservation of angular momentum in place 

of the conservation of tangential momentum for wind ac- 
cretion scenario, and considering the square and higher 
power terms of orbital eccentricity, the change equations 
of orbital elements have been recalculated. These equa- 
tions have been applied to quantitatively calculate the or- 
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bital elements of barium stars and to discuss the stellar 
wind accretion scenario. 

Our results show that during the wind accretion pro- 
cess, the system widens as it loses mass, the orbital semi- 
major axis A will increase, resulting in an increase in the 
orbital period, which is similar to the result of Boffin & 
Jorissen (1988), but the eccentricity is very nearly con- 
stant. In general, barium stars evolved from normal G, 
K giants, namely, barium stars formed by accreting the 
ejecta of the G, K giant companions of binary systems. 
If the accretion is through stellar wind, the orbital pe- 
riod will increase and eccentricity remains nearly constant. 
Thus, for the same eccentricity, the orbital periods of bar- 
ium stars are longer than those of normal G, K giants; 
while for the same orbital period, the maximum and aver- 
age eccentricity of barium stars are lower than the corre- 
sponding values of the normal G, K giants. Our results can 
explain quantitatively these relations. Namely, we can ex- 
plain the properties of (e,logP) diagrams of normal cluster 
giants and barium stars shown in Sect. 2. 

(1) Eq. (2) represents the change of orbital semi-major 
axis A, which shows that A will increase during the wind 
accretion process, resulting the increase of orbital periods. 
But the orbital eccentricity will be very nearly constant 
with the wind accretion (see Eq. (3) and Table 1). These 
results quantitatively explain the observational fact that 
the maximum eccentricities of barium stars is lower than 
the corresponding maximum eccentricities of normal gi- 
ants at any given orbital periods. 

(2) Regarding the max;imum eccentricity line of normal 
G, K giants as the initial condition of a binary system, 
the maximum eccentricity line of barium stars can be ob- 
tained through the wind accretion calculation, which can 
fit the observations. Moreover, the two maximum eccen- 
tricity lines are near parallel (see Fig. 1 and 2). 

(3) The average orbital eccentricity line of barium stars 
can be obtained through the wind accretion calculation, 
and can fit the observations. Namely, adopting rj = 0.685 
for Eq. (11) of BoSin et al. (1993) to obtain the average 
eccentricities of barium stars, the initial eccentricities of 
binary systems can be obtained using our wind accretion 
model, which nearly overlap the observed average values 
of normal G, K giants given by BoHin et al. (1993) (see 
Fig. 3). 

(4) For the low eccentricity systems, the observed max- 
imum orbital period of circular orbits is ^350 days for the 
cluster giants. Our result shows that this upper bound will 
lead to the corresponding maximum circular orbital period 
near 2000 days for barium stars using the wind accretion 
model, which can fit the observations. The threshold pe- 
riods of CH and S stars are very similar to that of barium 
stars. 

(5) Taking the highest orbital eccentricity 0.9 of nor- 
mal G, K giants as the initial condition, we can calculate a 

large eccentricity, up to 0.9, for barium stars. The reason 
may be that the orbital period increases, but the eccentric- 



ity is nearly constant during the wind accretion process, 
so that the long orbital period and high eccentricity bar- 
ium systems are formed. Our results can quantitatively 
explain the observed large eccentricity, 0.9, of some bar- 
ium systems. 

These fits between our results and the corresponding 

observations show that barium systems evolved from the 
normal G, K giants, and the distribution of orbital eccen- 
tricity with orbital period of barium stars can refiect the 
final orbital properties of binary systems through wind 
accretion for the binaries with longer period (P>1600 d). 

Moreover, this period value can be understood roughly 
by the results showed in Fig. 4. Fig. 4 shows the evolution 
of Roche radius as a function of mass ratio q = Mi/M^. 
The wind accretion scenario needs the larger Roche ra- 
dius than the stellar radii of the previous primary stars 
of the binary systems, AGB stars. If we also consider the 
radii of AGB stars: a typical AGB star with logTes = 3.5 
and Mboi = — 5 has a radius of 28Oi?0 (Jorissen & Boffin 
1992), the comparison between the evolution of the radii 
of the AGB stars and the evolution of the Roche radii will 
show that the Roche radii are larger than the correspond- 
ing stellar radii for the binary systems with longer final 
orbital period (P>1600 d) during the binary evolutions. 
This result reflect that the barium stars with longer or- 
bital period (P>1600 d) form through the wind accretion 
scenario. 

In addition, Liang et al. (2000) confirmed this orbital 
period value taking advantage of the observational heavy 
element abundances of barium stars. 

For the barium systems with shorter orbital periods, 
our wind accretion model is not sufficient to explain the 
observed properties. They may be formed through dynam- 
ical stability case C Roche lobe overflow, or common en- 
velope ejection, or some other scenarios (Jorissen et al. 
1995; Jorissen et al. 1998; Han et al. 1995 and references 
therein). Among the binary systems with shorter periods, 
the primary star may fill the Roche radius when it evolve 
to AGB phase, thus the mass-transfer process occur, but 
the previously strong wind mass loss has made the mass 
of primary star become very small so that mass ratio q 
(= M1/M2) is lower than a threshold value (/crit (^0.65, 
Pastetter & Ritter 1989), which causes the dynamically 
stable RLOF. While if the mass ratio q is greater than the 
(7crit when the primary AGB star fills its Roche lobe, the 
Roche lobe overflow will be dynamically unstable, leading 
to the formation of common envelope. Friction between 
the stars and the common envelope causes contraction of 
the binary orbit and release of energy. Most of the energy 
is used in the ejection of the common envelope. The binary 
after the ejection of common envelope comprises a white 
dwarf and a barium star with accreted neutron capture 
process elements— rich material. For such a binary, the or- 
bital period and eccentricity are lower. The barium star 
HD 77247 (P = 80.53 d, e = 0.09) can be explained by 
the common envelope scenario. 
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Appendix: derivation of Eqs. (2) and (3) 

For the binary system, the two components (an intrinsic 
AGB star with mass Mi, the present white dwarf star, and 
a main sequence star with mass M2, the present barium 
star) rotating around the mass core C, so the total angu- 
lar momentum is conservative in the mass core reference 
frame. If the mass exchange between the two components 
is wind accretion, the total angular momentum conserva- 
tion is shown by: 

A{iir^0) = LorliAMi + AM2) + r2v{AMi + AM2), (5) 

where fi is reduced mass, and r is the distance from M2 to 
Ml. ri and r2 are the distances from Mi, M2 to the mass 
core C respectively, lu {—2tt/P) is angular velocity, where 
P=2ttA'^{1 - e'^)i/h is the orbital period (Huang 1956). 
V is an additional effective velocity defined through the 
angular momentum variation in the direction of orbital 
motion of component 2. The first term on the right side 
of the equal-sign is the angular momentum lost by the 
escaping material and the second term is the additional 
angular momentum lost by the escaping material. 

For the binary system, according to Huang (1956), the 
changes of orbital elements, the orbital semi-major axis A 
and eccentricity e, are 
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According to the angular momentum conservation 
model, the A{r0) term can be obtained from the following 
equation: 

M1M2 . f M1M2 , 
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^ A{M-i_ + M2) ^^rG^M2A{Mi+ M2) 



M2 

then, wc can obtain 
reA{re) _ {rOf r 



E 



E 

Ar - 

E 



M-i{M-i + M2)^A2 

A(Mi + 



(9) 



Ar 



(10) 



M, 



AMi AM2 

Ml ^ M2 Ml + M2 

re A(Mi + M2) 

V 

E M2 



+ 



r^eG^M2A{Mi+M2) 
E Mi{Mi+M2)^Ai 

= 2(1 -e^) 



M2 



VoTh_ 



AMi AM2 A (Ml + M2) V 
Ml M2 
M2A(Mi + M2) 
Mi(Mi -hM2) 



2aA(Mi + M2) , MAr 
/{I e ) 



2(1 -e^)* 



Ml + M2 Er 
AMi AM2 A(i\/i + M2) V 



Ml 



M, 



Mo 



Vorb 



_ on 2a ^2A(Mi+M2) 

' Mi(Mi-fM2) 
-2(l-e^)^^y\+.y + 



2Ar 



Thus 



AE 



= 2(1 -e^)^ 



M1 + M2 A(l-e2)i 

A (Ml 



(11) 



AMi AM2 

+ 



M2 



Ml 



M, 



M2 



■forb 



_ gn 2a ^2A(Mi+M2) 

^ ^ Mi(Mi+M2) 

+ (2-2(l-e^)^)^(ML±^. 

^ ^ ^ M1 + M2 

The Ah/h term can be obtained from: 
A{^lr^e) = A{nh) 

^Mi + M2 ^ 

= Ai^^)h 



(12) 



^Mi + M2' 



M1M2 

mTTmz' 



■Ah 



and 
Ah 



wrlA{Mi + M2) + r2vA{Mi + M2) , (13) 



AMi AM2 ^ 2 + A(Mi + M2) v 



+ 



Ml M2 ' 2(1 -e2) 3 M2 
M2A(Mi + M2) A(Mi + M2) 



Mi(Mi+M2) 



Ml +M2 



Vorb 

(14) 



7 



Combining Eqs. (6), (7), (12) and (14), we can obtain 
the changes of orbital semi-major axis A and eccentricity 



AA 
~A 



-2(1 



e2)5 



AMi _^ AM2 



AMi + AM2 V 



Ml 



M2 



Mo 



Vorh 



, na ^2(AMi+AM2) 
^ ^ Mi(Mi+M2) 
,,AMi + AM2 



+ [2(1 -e^)^ 



Ml + i\f2 



(15) 



sAe 



1-e^ 



= [l-(l-e2)^] 



AMi AAf2 AMi + AM2 



Ml 



M2 



Ml + M2 



2 1 M2(AMi + AM2) 
^ ^ ^ ^ ^ Mi(Mi +M2) 
3e2 AMi + AM2 

2(l_e2)i M2 l^orb' 



(16) 



For obtaining the above equations, averaging over the 
orbital period then amounts to using the following mean 
values: 

_ 2tt 
^ r PA{1 - e2) 



<r>= A{1 + -e^) 



< r9 >= 



<0>-- 



27rA(l-e2)i 



P 



27r 
"P 



P 



1 1 

< - >= -r 
r A 
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A2(l-e2)i 



